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ABSTRACT 

Context. Irregular galaxies are usually smaller and less massive than their spiral, SO, and elliptical counterparts. Radio observations 
indicate that a magnetic field is present in irregular galaxies whose value is similar to that in spiral galaxies. However, the conditions 
in the interstellar medium of an irregular galaxy are unfavorable for amplification of the magnetic field because of the slow rotation 
and low shearing rate. 

Aims. We investigate the cosmic-ray driven dynamo in the interstellar medium of an irregular galaxy. We study its efficiency under 
the conditions of slow rotation and weak shear. The star formation is also taken into account in our model and is parametrized by the 
frequency of explosions and modulations of activity. 

Methods. The numerical model includes a magnetohydrodynamical dynamo driven by cosmic rays that is injected into the interstellar 
medium by randomly exploding supernovae. In the model, we also include essential elements such as vertical gravity of the disk, 
differential rotation approximated by the shearing box, and resistivity leading to magnetic reconnection. 

Results. We find that even slow galactic rotation with a low shearing rate amplifies the magnetic field, and that rapid rotation with 
a low value of the shear enhances the efficiency of the dynamo. Our simulations have shown that a high amount of magnetic energy 
leaves the simulation box becoming an efficient source of intergalactic magnetic fields. 

Key words. MHD - ISM: magnetic fields - Galaxies: irregular - Methods: numerical 



1. Introduction 

Irregular galaxies have lower masses than typical spirals and el- 
lipticals. In addition, they have irregular distributions of the star- 
forming regions, and rotations that are slower than spiral galax- 
ies by half an order of magnitude (Gallagher & Hunter 1984). 
The rotation curves of irregular galaxies are non-uniform and 
have a weak shear. 

Radio observations of magnetic fields in spiral galaxies in- 
dicate that their magnetic fields have strong ordered (1-5 fiG) 
and random (9-15 yuG) components (Beck 2005). A plausible 
process fueling the growth in the magnetic energy and flux 
of these galaxies is magnetohydrodynamical dynamo (Widrow 
2002; Gressel et al. 2008). The vital conditions required for the 
dynamo to effectively amplify the magnetic field are rapid rota- 
tion and shear. In irregular galaxies, both quantities seem to be 
too low to initiate efficient dynamo action. In contrast, the ob- 
servations of magnetic field in irregular galaxies indicate that 
these galaxies could have strong and ordered magnetic fields 
(e.g., Chyzy et al. 2000, 2003; Kepley et al. 2007; Lisenfeld et 
al. 2004). 

The most spectacular radio observations to date of irregu- 
lars were those performed for the galaxy NGC 4449 (Chyzy et 
al. 2000). The total strength of its magnetic field is about 14 fiG 
with a ordered component reaching locally values of 8 fiG. These 
are similar to the intensities observed for large spirals. A high 
number of H n regions and slow rotation is also observed with 
quite large velocity shear (Valdez-Gutierrez et al. 2002). The ra- 
dio observations of H i around the galaxy indicate that this object 
is embedded in two large H i systems that counter-rotate with re- 



spect to the optical part of this galaxy (Bajaja et al. 1994; Hunter 
et al. 1998, 1999). In addition to these Hi clouds, NGC 4449 
contains an unusual ring of H i in the outer part of the optical 
disk (Hunter et al. 1999). This complicated topology of the Hi 
velocity field could help in achieving efficient magnetic field am- 
plification (see Otmianowska-Mazur et al. 2000). 

Chyzy et al. (2003) found that two other irregular galaxies, 
NGC 6822 and IC 10 are also magnetized. The former has a very 
low total magnetic field weaker than 5 /iG, a small number of 
Hn regions, and almost rigid rotation (see Sect. [2}. These prop- 
erties are directly related to the efficiency of the dynamo pro- 
cess in galaxies (see Otmianowska-Mazur et al. 2000; Hanasz et 
al. (2006, 2009) that weakly amplifies the magnetic field in this 
galaxy. The irregular IC 10 has a total magnetic field strength 
that varies between 5 and 15 fiG with no ordered component. 
Observations performed by Chyzy et al. (2003) indicate that the 
total magnetic field is correlated with the number of H n regions. 
The number of the regions is higher than in NGC 6822, and the 
rotation of IC 10 has a partly differential character (see Sect. [2j. 
Both conditions lead to more rapid magnetic field amplification 
than for NGC 6822. As for NGC 4449, IC 10 is embedded in 
a large cloud of Hi, which counter-rotates with respect to the 
inner disk (Wilcots & Miller 1998). 

Klein et al. (1993) inferred that the Large Magellanic Cloud 
(LMC) has a large-scale magnetic field that has the shape of 
a trailing spiral structure, similar to normal spiral galaxies. It is 
possible that the amplification of the magnetic field is connected 
to the differential rotation of this galaxy present beyond a cer- 
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Table 1. Main properties of NGC 4449, NGC 6822, and IC 10 





NGC 4449 


NGC 6822 


IC 10 


Type 


IBm 


IB(s)m 


dlrr 


Distance [Mpc] 


4.21 


0.50 


0.66 


Diameter [kpc] 


5.73 


1.71 


1.27 


SFR 


high 


low 


high 


<B to t>LuG] 


5-15 


< 5 


5-15 


<B reg > [juG] 


=! 8 


< 3 


< 3 



In the following rows we present: the morphological type of the object 
(LEDA), the distance (Karachentsev et al. 2004), the size calculated 
from £>25 (LEDA) and distance, the star formation rate (SFR), and the 
results of the radio analysis of Chyzy et al. (2000, 2003). 



tain radius (Klein et al. 1993; Luks & Rohlfs 1992; Gaensler et 
al. 2005). 

We note that polarized radio emission is detected in the ir- 
regular galaxy NGC 1569. This galaxy has a very high star for- 
mation rate (Martin 1998) and exhibits bursts of activity in its 
past (Vallenari & Bomans 1996). The radio observation of this 
galaxy by Lisenfeld et al. (2004) found that the galaxy has large- 
scale magnetic fields in the disk and halo. Furthermore, they 
found that their data agree that a convective wind could allow 
for escape of cosmic-ray electrons in to the halo. These observa- 
tions are the main reason for undertaking our CR-driven dynamo 
calculations in irregular galaxies. In addition to Lisenfeld et al. 
(2004), the radio observations of Kepley et al. (2007) showed 
that the large-scale magnetic arms visible in NGC 1569 are 
aligned perpendicularly to the disk and that the northern part of 
the disk of the galaxy is inclined at a different angle. 

Kronberg et al. (1999) realized that dwarf galaxies (apart 
from their low masses) could serve as an efficient source of gas 
and magnetic fields in the intergalactic medium (IGM) during 
their initial bursts of star formation in the early Universe. In the 
case of star-forming dwarf galaxies, we expect that the dominant 
driver of a galactic wind are cosmic rays, in contrast to large spi- 
rals, such as the Milky Way, where the thermal driving is most 
significant (Everett et al. 2008). Therefore, we applied the model 
of the CR-driven dynamo to the interstellar medium and con- 
ditions of an irregular galaxy and try to find how much of the 
magnetic energy can be expelled from the dwarf galaxies to the 
IGM. 

Many questions about magnetic field amplification in irreg- 
ular galaxies remain unresolved. The physical explanation of 
this process is difficult to establish because these galaxies rotate 
slowly, almost like a solid body. In this paper, we check how our 
model of cosmic -ray driven dynamo, which effectively describes 
spiral galaxies (Hanasz et al. 2004, 2006, 2009), can be applied 
to irregular galaxies. In the present numerical experiment we at- 
tempt to answer how the model input parameters observed in ir- 
regulars (small gravitational potential, gas density, low rotation, 
and small shear) influence the magnetic fields within them. We 
have not taken into account the magnetic field possibly injected 
by stars. We plan to study this in the future. We found that in 
certain conditions achievable for irregulars it is possible to have 
efficient magnetic field amplification. 



2. Observations of irregular galaxies 

To study properties of the irregulars and determine the input pa- 
rameters for our simulations, we use observations of NGC 4449, 



NGC 6822, and IC 10 acquired by Chyzy et al. (2000, 2003). 
The main properties of these objects are presented in Table Q] 

From the rotation curves of IC 10 (Fig. Q] top panel, solid 
line) obtained by Wilcots & Miller 1998, NGC 6822 (Fig. [Q 
top panel, dashed line) obtained by Weldrake et al. 2003, and 
NGC 4449 (Fig. Q] top panel, dot-dashed line) by Valdez- 
Gutierrez et al. 2002, we computed the angular velocity and 
shearing rate of each galaxy (see Sect. 15.1b . For the two first 
galaxies, we use H i data. In the case of NGC 4449, we restricted 
our analysis to the internal region and used Ha data, because of 
its very complex velocity pattern. 

In the velocity pattern of the IC10, we can see a central part 
with a solid-body rotation, which flattens to a constant value 
v ro t - 30 km s _1 at r — 352 pc. The NGC 6822 rotation curve is 
a monotonically increasing function with a square root slope and 
the highest value v rot - 60 km s _1 at r — 5.7 kpc. The rotation 
curve of NGC 4449 is highly disturbed and reaches a maximum 
value of Vrot - 40 km s _1 at radius of 2 kpc. 



3. Description of the model 

The CR-driven dynamo model consists of the following ele- 
ments (based on Hanasz et al. 2004, 2006): 

(1) The cosmic ray component is a relativistic gas described by 
a diffusion-advection transport equation. Typical values of 
the diffusion coefficient are (3 -f 5) x 10 28 cm 2 s _1 (see Strong 
et al. 2007) at energies of around 1 GeV, although in our 
simulations we use reduced values (see Sect. 14. 11 1. 

(2) Anisotropic diffusion of CR. Following Giacalone & Jokipii 
(1999) and Jokipii (1999), we assume that the CR gas dif- 
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Fig. 1. Observational rotation characteristics of IC 10, 
NGC 6822, and NGC 4449. We present, from top to bot- 
tom: the rotation curves (references in Sect. 0, the calculated 
angular velocity, and the computed shear parameter q (for 
details see Sect. 15. U respectively. The shaded region marks the 
range of parameters presented in this paper. 
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fuses anisotropically along magnetic field lines. The ratio of 
the perpendicular to parallel CR diffusion coefficients sug- 
gested by the authors is 5%. 

(3) Localized sources of CR. In the model, we apply the random 
explosions of supernovae in the disk volume. Each explo- 
sion is a localized source of cosmic rays. The cosmic ray 
input of individual SN remnant is 10% of the canonical ki- 
netic energy output (10 51 erg) and distributed over several 
subsequent time steps. 

(4) Resistivity of the ISM to enable the dissipation of the small- 
scale magnetic fields (see Hanasz et al. 2002 and Hanasz & 
Lesch 2003). In the model, we apply the uniform resistivity 
and neglect the Ohmic heating of gas by the resistive dissi- 
pation of magnetic fields. 

(5) Shearing boundary conditions and tidal forces following the 
prescription by Hawley, Gammie & Balbus (1995), are im- 
plemented to reproduce the differentially rotating disk in the 
local approximation. 

(6) Realistic vertical disk gravity following the model by 
Ferriere ( 1 998) modified by reducing the contribution of disk 
and halo masses by one order of the magnitude, to adjust the 
irregular galaxy environment. 



We apply the following set of resistive MHD equations: 



dp 

| + V.(pV) = 0, 

^ + V • (e\) = -p(V ■ V), 

d\ 1 / B 2 



+ 



B VB 

4np 



■ 2Q x V + 2qQ. 2 xe x 



dB 

— = Vx(VxB) + i/aB, 
ot 

p = (y- y = 5/3, 



(1) 

(2) 

(3) 

(4) 
(5) 



where q = -din Q./dlnR is the shearing rate, R is a galactocen- 
tric radius, 77 represents the ISM resistivity, y is the adiabatic in- 
dex of thermal gas, p a is the cosmic-ray pressure, and the other 
symbols have their usual meaning. In the equation of motion, the 
term Vp cl is included (see Berezinskii et al. 1990). The thermal 
gas is approximated by an adiabatic medium. 

The cosmic ray component is an additional fluid described 
by the diffusion-advection equation (see e.g., Schlickeiser & 
Lerche 1985) 



de a 
dt 



+ V(e CI V) = V(KVe c[ ) - p CT (V -V) + Q : 



SN, 



where Qsn is the source term of the cosmic-ray energy den- 
sity injected locally from randomly exploding SN remnants. The 
cosmic-ray fluid is described by an adiabatic equation of state 
with adiabatic index y cr : 

Per = (Tcr - lKr, Tor = 14/9. (7) 

The K is an diffusion tensor described by the formula: 

K ij = K x 8 i j + (K l[ -K x )n i n j , n^Bi/B, (8) 

adopted following the argumentation of Ryu et al. (2003). 

The vertical gravitational acceleration is taken from Ferriere 
(1998). We reduced both contributions of disk and halo by a fac- 
tor of 10, the scale length of the exponential disk to La — 2 kpc, 
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Fig. 2. Example slices of a domain taken from simulation 
R.OIQI at t — 660 Myr. On the slices, the Parker loop is pro- 
duced by cosmic rays from supernovae explosions. 



and scale length of halo to L# = 1 kpc. In our computations, we 
incorporated the formula : 



R + L l 2 
? Z (R,Z) = (1.7- 10~ 10 cms~ 2 ) ' 



+ (4.4. l(T 10 cms" 2 )exp 



R + L 2 H \ 1 kpc 
R-RA Z 



Ld ) y/Z 2 + (0.2 kpc) 2 



, (9) 



where R* is the distance of the origin of the simulation box from 
the galactic center and Z is the height above the galactic mid 
plane. 



(6) 4. Model setup and parameters 



4.1. Model setup 

The 3D cartesian domain size is 0.5 kpc x 1 kpc x 8 kpc in x, y, z 
coordinates corresponding to the radial, azimuthal, and vertical 
directions, respectively. The grid size is 20 pc in each direction. 
The boundary conditions are sheared-periodic in x, periodic in 
y, and an outflow in z direction. The domain is placed at the 
galactocentric radius R* = 2 kpc. In Fig. [2] we present example 
slices through the simulation domain. The left panel shows the 
CR energy density with the magnetic field vectors and the right 
panel shows the gas density with velocity vectors. 

The positions of SNe are chosen randomly with a uniform 
distribution in the xy plane and a Gaussian distribution in the ver- 
tical direction. The scaleheight of SN explosions in the vertical 
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Table 2. List of models 



Model 


Ll 


1 


J 




T„ 




[Myr -1 ] 




[kpc -2 Myr -1 ] 


[Myr] 


[Myr] 


R.OIQI 


0.01 


1 


10 


200 


20 


R.02Q1 


0.02 


1 


10 


200 


20 


R.03Q1* 


0.03 


1 


10 


200 


20 


R.04Q1 


0.04 


1 


10 


200 


20 


R.05Q1 1 


0.05 


1 


10 


200 


20 


R.01Q0 


0.01 





10 


200 


20 


R.01Q.5 


0.01 


0.5 


10 


200 


20 


f\ 1 /" "\ 1 ft 

R.OIQI 


0.01 


1 


10 


200 


20 


R.01R1.5 


0.01 


1.5 


10 


200 


20 


R.05Q0 


0.05 





10 


200 


20 


R.05Q.5 


0.05 


0.5 


10 


200 


20 




n n*> 

U.Uj 


i 
i 


i n 
l u 


inn 
ZUU 


in 
zu 


r> d 1 c 
K.UjKI.j 


n n<* 
U.U3 


i.j 


i n 
1U 


inn 
ZUU 


ZU 


CE1D ATA C 


U.UJ 


U.J 


J 


inn 
ZUU 


on 
ZU 


Con ( \~) f \ 1 


0.03 


1 


3 


200 


20 


SF10R.03Q.5 


0.03 


0.5 


10 


200 


20 


SF10R.03Q1 


0.03 


1 


10 


200 


20 




n n^i 

U.UJ 


n 

U.J 


^n 

JU 


inn 

ZUU 


on 

ZU 


SF30R.03Q1 


0.03 


1 


30 


200 


20 


Ml 0/1 00 


0.03 


0.5 


10 


100 


10 


M20/200 f/ 


0.03 


0.5 


10 


200 


20 


M50/100 


0.03 


0.5 


10 


100 


50 


Ml 00/200 


0.03 


0.5 


10 


100 


200 


M100/100 


0.03 


0.5 


10 


100 


100 


FIRST 


0.03 


0.5 


2.5 


2000 


50 



Subsequent columns show: the model name, the angular velocity D, 
the shearing parameter q, the frequency of SN explosions /, the period 
of SNe modulation T p , a nd the duration of SNe activity in one the 
period T a . See Sect. l4.2l for details. The horizontal lines distinguish 
between different simulation series. Models with the same superscript 
( a , *, c , and d ) point to the same experiments, but are written for clarity. 



direction is 100 pc, and the CR energy that originates in an ex- 
plosion is injected instantaneously into the ISM with a Gaussian 
radial profile (rsN = 50 pc). In addition, the SNe activity is mod- 
ulated during the simulation time by a period T p and an activity 
time T a . 

The applied value of the perpendicular CR diffusion coeffi- 
cient is K ± = 10 3 pc 2 Myr" 1 = 3 x 10 26 cm 2 s -1 and the parallel 



27 



The diffusion 



one is K\\ = 10 4 pc 2 Myr' = 3 x 10 
coefficients are 10% of realistic values because of the simulation 
timestep, which becomes prohibitively short when the diffusion 
is too high. 

The initial state of the system represents the magnetohydro- 
static equilibrium with the horizontal, purely azimuthal mag- 
netic field with Pmag/Pgas = 10 -4 , which corresponds to the 
mean value of magnetic field in the simulation box of 5 nG. 
Magnetic diffusivity 77 is set to be 100 pc 2 Myr -1 , which corre- 
sponds to 3x 10 25 cm 2 s -1 in cgs units (Lesch 1993). The column 
density of gas is g gas = 6 x 10 20 cm -2 (taken from observations, 
see Gallagher & Hunter 1984) and the initial value of the isother- 
mal sound speed is set to be c, JO = 7 km s -1 . 



4.2. Model parameters 

We present the results of four simulation series corresponding to 
different sets of the CR-dynamo parameters. Details of all com- 
puted models are shown in Table [2] The model name consists 
of a combination of four letters: R, Q, SF and M followed by 
a number. The letter R means the angular velocity (rotation), Q is 



the shearing rate, SF is the supernova explosion frequency and 
M represents for its modulation during the simulation time, and 
the numbers determine the value of the corresponding quantity. 
Only the modulation symbol is followed by two numbers, the 
first corresponding to the time of the SNe activity and the sec- 
ond to a period of modulation. Values of the parameters are given 
in the following units: angular velocity in Myr -1 , supernova ex- 
plosion frequency in kpc -2 Myr -1 , and the modulation times in 
Myr. For example, a model named "R.OIQI" denotes a simula- 
tion where Q. = 0.01 Myr -1 and q = 1, and the name "M50/100" 
denotes an experiment with T a = 50 Myr and T p = 100 Myr. The 
last model in Table [2] named FIRST, points to an experiment, in 
which only during the first 50 Myr supernovae are active and 
after that time CR injection stops. 

5. Results 

5.1. Shear parameter q obtained from observations 

The shearing rate parameter q (defined in Sect. is calculated 
numerically from the observational rotation curves using a sec- 
ond order method 



1 Ri 
qt = l . . 

Vi \Ri+\ 



v i+ i-Vi-i 1 Vj+i - 2v t + v,-_i 
h 

Ri-l 2 Ri+i - 



(10) 



applied to the radial velocities v, measured at Rj of the observed 
rotation curve. Calculations are performed only where the rota- 
tion curve is smooth enough, because of the enormous velocity 
fluctuations and low spatial resolution, which cause large dis- 
persions in our results. The estimated shearing parameters from 
observational rotation curves are presented in Fig. Q] Different 
values of the parameter q, correspond to the following interpre- 
tations: when q < 0, the rotation velocity increases faster than 
a solid body; when q — 0, we have solid body rotation; when 
< q < 1, the rotation velocity increases slower than a solid 
body; q = 1 relates to a flat rotation curve; for q > 1, the az- 
imuthal rotation decreases with R. We found that the shearing 
rates are high in all three galaxies and due to the large variations 
in the rotation curves, q changes rapidly. However, in the case 
of NGC 6822, q gradually increases from to 1.5 with galac- 
tocentric distance. For the galaxies IC 10 and NGC 4449, the 
calculated local shearing rates vary from -1.5 to 3. This scatter 
in the results is caused by the large fluctuations in the measured 
rotation velocities. 



5.2. The magnetic field evolution 

We study dependence of the magnetic field amplification on the 
parameters describing the rotation curve, namely, the shearing 
rate q and the angular velocity O. The evolution in the total 
magnetic field energy Eb and total azimuthal flux for dif- 
ferent values of Q is shown in Fig. [3] left and right panel, re- 
spectively. Models with higher angular velocities, starting from 
0.03 Myr -1 , initially exhibit exponential Eb growth and after 
about 1 200 Myr, the process saturates (see Sect.|6]for the discus- 
sion). The saturation values of magnetic energy for these three 
models are similar and Eb exceeds the value 10 4 in the nor- 
malized units. The magnetic energy in the models R.OIQI and 
R.02Q1 grows exponentially during the whole simulation and 
does not reach the saturation level. The final Eb for R.02Q1 is 
around 4 x 10 3 and for the slowest rotation (R.OIQI) in our 
sample is only 20. The total azimuthal magnetic flux evolu- 
tion (Fig. [3] right) shows that a higher angular velocity pro- 
duces a higher amplification. The azimuthal flux for models with 
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Model 


log£™" 


log E°<>< 


JiOiit j fend 
B ' B 


(B) 


t e 


T n 


Description 










\»G] 


[Myr] 


[Myr] 




R.01Q1" 1 


1.37 


0.07 


0.05 


0.068 


1219 


628 


slow rotation 


R.02Q1 


3.55 


2.02 


0.03 


0.835 


440 


314 


slow rotation 


R.03Q1 6 


4.03 


2.89 


0.07 


1.206 


375 


209 


medium rotation 


R.04Q1 


4.17 


3.29 


0.13 


1.285 


346 


157 


fast rotation 


R.05Q1 1 


3.95 


3.43 


0.30 


1.120 


509 


125 


fast rotation 


R.01Q0 


-1.73 


-0.42 


20.48 


0.001 


- 


628 


low shear 


R.01Q.5 


0.06 


-0.35 


0.39 


0.011 


- 


628 


medium shear 


R.01Q1" 


1.37 


0.07 


0.05 


0.068 


1219 


628 


medium shear 


R.01Q1.5 


1.26 


0.15 


0.08 


0.056 


1232 


628 


high shear 


R.05Q0 


-1.00 


0.04 


10.78 


0.003 


- 


125 


low shear 


R.05Q.5 


4.03 


3.38 


0.22 


1.168 


363 


125 


medium shear 


R.05Q1 1 


3.95 


3.43 


0.30 


1.120 


509 


125 


medium shear 


R.05Q1.5 


3.81 


3.07 


0.18 


0.934 


375 


125 


high shear 


SF3R.03Q.5 


3.14 


2.23 


0.12 


0.190 


506 


209 


low SFR 


SF3R.03Q1 


3.28 


2.28 


0.10 


0.243 


444 


209 


low SFR 


SF10R.03Q.5'' 


3.87 


2.28 


0.03 


0.991 


403 


209 


medium SFR 


SF10R.03Q1* 


4.03 


2.89 


0.07 


1.206 


375 


209 


medium SFR 


SF30R.03Q.5 


2.79 


2.77 


0.96 


0.190 


629 


209 


high SFR 


SF30R.03Q1 


3.30 


3.05 


0.57 


0.243 


547 


209 


high SFR 


M10/100 


3.78 


2.73 


0.09 


0.833 


422 


209 




M20/200 rf 


3.87 


2.28 


0.03 


0.243 


403 


209 




M50/100 


4.05 


2.92 


0.07 


1.001 


404 


209 




Ml 00/200 


4.11 


2.67 


0.04 


1.352 


393 


209 




M100/100 


3.86 


2.72 


0.07 


0.751 


422 


209 


constantly exploding SNe 


FIRST 


2.37 


1.33 


0.09 


0.113 


572 


209 


SN activity during first 50Myr 



The subsequent columns show: the model name, the mean value of total magnetic energy Eg" 1 over past 50 Myr, the total outflow of magnetic 
energy during whole simulation (see Eq. ll2b . the ratio of these two quantities, the final mean value of the magnetic field in a disc midplane (6) 
(|z| < 20 pc), the e-folding time of magnetic flux increase t e , galaxy revolution timescale T n , and a short description of a model. Superscripts are 
explained in Table[2] Values of magnetic field energy are normalized to the initial value. 



O > 0.02 Myr 1 exceeds the value 10 2 . Model R.01Q1 does not 
enhance the azimuthal flux at all. 

In Fig. [4] we present results for models with different shear- 
ing rate values. The evolution of Eb and B^ in models R.05Q.5 
and R.05Q1.5 follows the evolution of model R.05Q1, which is 
described in the previous paragraph. Similar behavior is noted 
for R.01Q1.5 and R.01Q1, but the model R.01Q.5 alone sustains 
its initial magnetic field. In the case of models with no shear 
(R.01Q0, R.05Q0), the initial magnetic field decays. 

We check how the frequency and modulation of SNe influ- 
ence the amplification of magnetic fields. The evolution in total 
magnetic field energy and total azimuthal flux for different su- 
pernova explosion frequencies are shown in Fig. [5] left and right 
respectively. The total magnetic energy evolution for all mod- 
els is similar, but in the case of the azimuthal flux we observe 
differences between the models. The most efficient amplifica- 
tion of fi appears for SF10R.03Q.5 and SF10R.03Q1, and for 
other models the process is less efficient. In addition, for mod- 
els SF30R.03Q.5 and SF30R.03Q1, we observe a turnover in 
magnetic field direction. The results suggest that the dynamo 
requires higher frequencies of supernova explosions to create 
more regular fields, although, if the explosions occur too fre- 
quently, this process is suppressed because of the overlapping 
turbulence. The analysis of the M models (Fig. [6]l shows that 
the dynamo process depends on the duration of the phase when 
supernova activity switches off. The fastest growth of magnetic 
field amplification occurs for models M100/200 and M50/100 in 
which periods of SN activity occupy half of the total modulation 
period. The amplification is apparently weaker in cases of short 
SN activity periods (Ml 0/1 00, M20/200) and continuous activ- 



ity (M100/100). In all M models, the final E B reaches a value 
of the order of 10 4 . For evolution, we found that the mag- 
netic flux in the model M50/100 increases exponentially and sat- 
urates after 1 300 Myr. Similar behavior is exhibited by the mod- 
els M10/100 and M100/100 but the saturation times occur after 

1 700 Myr and the growth is slower than in the previous case. 
The model Ml 00/200 after exponential growth at t = 1 650 Myr 
probably begins to saturate, but to quantify this exactly the simu- 
lation should continue. The model M20/200 grows exponentially 
and does not appear to saturate. 

In the case of the model FIRST (Fig. [6]), we found that af- 
ter about 8 galaxy revolutions the growth in Eb and B^ stops. 
The total magnetic field energy increase exponentially and af- 
ter reaching a maximum at t — 1 400 Myr, it exceeds the value 

2 x 10 2 , whereas the azimuthal flux saturates after 1 600 Myr and 
afterwards starts to decay gently. 



5.3. Magnetic field outflow 

To measure the total production rate of magnetic field energy 
during the simulation time, we calculate the outflowing E° B l " 
through the xy top and bottom domain boundaries. To estimate 
the magnetic energy loss, we compute the vertical component of 
the Poynting vector 



S z = (B x v x + B yVy )B z - (B 2 X + B 2 y )v z . (11) 



6 



H. Siejkowski et al.: CR driven dynamo in the ISM of irregular galaxies 




R.01Q0 

R.01Q.5 

R.01Q1 

R.Q1Q1.5 

R.05Q0 

R.05Q.5 

R.05Q1 

R.05Q1.5 




.1 100 



20x R.01Q0 
20X R.01Q.5 
20x R.01Q1 
20x R.01Q1.5 
R.05Q0 
R.05Q.5 
R.05Q1 
R.0501 .5 



1000 
le [Myr 



1000 
le [Myr] 



Fig. 4. Evolution of the total magnetic energy Eb (left panel) and the total azimuthal flux (right) for models with different 
shearing rate and rotation. Both quantities are normalized to the initial value. For models with Q = 0.01 Myr 1 , the values has 
been multiply by factor 20 - proper y-axis for these plots is on the right side of the frame. 



This value is computed in every cell belonging to the top and 
bottom boundary planes and then integrated over the entire area 
and time: 



At. 



(12) 



models with q — 0, in which the dynamo does not operate, the 
outflowing energy originate only from the initial condition for 
the magnetic field. The results show that the outflowing mag- 
netic energy is substantial (see TableO suggesting, that irregular 
galaxies because of their weaker gravity can be efficient sources 
of intergalactic magnetic fields. 



where k mm and k max refer to the bottom and top boundary respec- 
tively, A z is the vertical dimension of single cell, t is the simula- 
tion time, and At is the timestep. For models with a low dynamo 
efficiency most of the initial magnetic field energy is transported 
out of the simulation box. In some cases (i.e., all models ex- 
cept R.01Q0 and R.05Q0), we find that the energy loss E° B ' is 
comparable to the energy remaining inside the domain E B ' d . In 
these models, the ratio E B " /E e B ' d varies from 0.03 to 0.96 and 
is highly dependent on the supernova explosion frequency. For 



6. Discussion 

The most effective magnetic field amplification that we have 
found is that in the model R.04Q1, which we associate with the 
galaxy NGC 4449. This galaxy has the highest star formation 
rate in our sample of three irregulars. The rotation is rapid, reach- 
ing 40km/s, and, for a wide range of radii, the shear is strong. 
The numerical model predicts an effective magnetic field am- 
plification and NGC 4449 indeed hosts the strongest magnetic 
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Fig. 5. Evolution of the total magnetic energy Eg (left panel) and the total azimuthal flux B^, (right) for models with different 
supernova explosion frequency and shearing rate. Both quantities are normalized to the initial value. 
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Fig. 6. Evolution of the total magnetic energy Eg (left panel) and the total azimuthal flux (right) for models with different times 
of supernovae modulation. Both quantities are normalized to the initial value. 



field among the irregulars, both in terms of its total and ordered 
component of 14 /iG and 8 /iG, respectively (Chyzy et al. 2000). 

The next galaxy IC 10 forms stars at a lower rate than 
NGC 4449. The shear is strong and it is a rapid rotator. We can 
compare this galaxy to our model R.05Q1, where we see the 
fastest initial growth of the total magnetic energy, but the final 
value is smaller than that in the case of slower rotation. The to- 
tal azimuthal flux evolves in a complex way with a reversal in 
the mean magnetic field direction. This may indicate that be- 
cause of its relatively rapid rotation and small size, instabilities 
can evolve faster. Separate instability domains can mix (over- 
lap) with each other resulting in a chaotic though still amplified 
magnetic field. Consequently IC 10 exhibits a strong total mag- 
netic field of 5-15//G (as estimated by Chyzy et al. 2003). We 
notice that by increasing the rotation speed, the amount of mag- 
netic energy expelled from the galaxy grows (see Table[3]l. IC 10 



has a relatively low mass and its shallow gravitational potential 
makes the escape of its magnetized ISM easier. 

NGC 6822 forms stars at the slowest rate in our sample. It is 
also the slowest rotator. The rotation is almost rigid in its cen- 
tral part (out to ~0.5kpc) gradual becoming differential at larger 
galactocentric distances but the calculated shearing rate remains 
small. We can explain its weak magnetic field of lower than 5 /iG 
(Chyzy et al. 2003) by comparing with our model FIRST: a sin- 
gle burst of star-forming activity in the past followed by a long 
(lasting until present) period of almost no star-forming activity. 
In this model, the magnetic field, amplified initially, fades since 
the star formation stops. This star-forming activity was analyzed 
for spiral galaxies by Hanasz et al. (2006), who measured a lin- 
ear growth in the magnetic field. We can explain this by using 
a shorter simulation time (by about a factor of two) than in our 
case, but it may indicate that in irregulars the magnetic field is 
more easily expelled from the galaxy. 
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Our models, for which we measure amplification in Eg and 
Z?^ during the simulation, produce a mean magnetic field of order 
1-0.5 fiG (Table [3} within a disc volume. Models with slower 
growth of magnetic field reach values of (B) around tens of nG, 
and models with no dynamo action diffuse the initial magnetic 
field outside the simulation box. 

In Table[3] we present the average e-folding time of the mag- 
netic flux increase t e and the galactic revolution period Tci- The t e 
of most models is in between 300 and 600 Myr. For spiral galax- 
ies, Hanasz et al. (2006, 2009) found that the e-folding timescale 
is about 150-190 Myr. The difference between spirals and irreg- 
ulars is probably caused by rotation, which is much more rapid 
in spirals. 

In most of our models, large fractions of the magnetic field 
are expelled out of the computational domain - almost 20%- 
30% of the magnetic energy maintained in the galaxy. In gen- 
eral more rapid rotation and a high SNe rate make it easier for 
the magnetized medium to escape. However, for higher shear 
rate, the share of the expelled magnetic field is lower. The 
optimal set of parameters, from this point of view, is repre- 
sented by the model R.05Q1, which we relate to IC 10. In the 
other two galaxies, the expelled field is also high - about 10%. 
Models with excessive star formation increase this fraction to 
60% (SF30R.03Q.5) or even 96% (SF30R.03Q1). Therefore, the 
irregular galaxies, in particular compact and intensively forming 
stars such as IC 10, are an important source of magnetic field in 
intergalactic and intracluster media, as predicted by Kronberg et 
al. (1999). 

For most of our models we found that the value of the mag- 
netic field strength in the vicinity of a galaxy (at z = 4 kpc) 
is about 30-200 nG. Only models with low magnetic-field pro- 
duction rates produce negligible magnetic fields at this height. 
This area is the highest point in our simulation domain above 
the galactic midplane and can be considered as a transition re- 
gion between the ISM and the IGM. Hence, the magnetic field 
strengths in the models can be an upper limit to the values in the 
IGM region. Our estimates are in an agreement with previous 
studies, including Ryu et al. (1998), who demonstrated that in 
largescale filaments, magnetic fields of about 1 fiG may exists, 
Kronberg et al. (1999), who calculated that on Mpc scales the av- 
erage magnetic field strength is about 5 nG, and Gopal-Krishna 
& Wiita (2001), who showed that radio galaxies can seed the 
IGM with a magnetic field of the order 10 nG during the quasar 
era. However, to obtain realistic profile or even the maximum 
possible range of expelled magnetic field in the case of dwarf 
galaxies we should take into account the interaction between the 
IGM and ISM (pressure), which is not included in our model. 
We plan to extend our research in this respect in future work. 



7. Conclusions 

We have described the evolution in the magnetic fields of ir- 
regular galaxies in terms of a cosmic-ray driven dynamo. Our 
cosmic -ray driven dynamo model consists of (1) randomly ex- 
ploding supernovae that supply the CR density energy, (2) shear- 
ing motions due to differential rotation, and (3) ISM resistivity. 
We have studied the amplification of magnetic fields under dif- 
ferent conditions characterized by the rotation curve (the angu- 
lar velocity and the shear) and the supernovae activity (its fre- 
quency and modulation) typical of irregular galaxies. We have 
found that: 



- in the presence of slow rotation and weak shear in irregular 
galaxies, the amplification of the total magnetic field energy 
is still possible; 

- shear is necessary for magnetic field amplification, but the 
amplification itself depends weakly on the shearing rate; 

- higher angular velocity enables a higher efficiency in the CR- 
driven dynamo process; 

- the efficiency of the dynamo process increases with SNe ac- 
tivity, but excessive SNe activity reduces the amplification; 

- a shorter period of halted SNe activity leads to faster growth 
and an earlier saturation time in the evolution of azimuthal 
magnetic flux; 

- for high SNe activity and rapid rotation, the azimuthal flux 
reverses its direction because of turbulence overlapping; 

- because of the shallow gravitation potential of an irregular 
galaxy, the outflow of magnetic field from the disk is high, 
suggesting that they may magnetize the intergalactic medium 
as predicted by Kronberg et al. (1999) and Bertone et al. 
(2006). 

The performed simulations indicate that the CR-driven dy- 
namo can explain the observed magnetic fields in irregular 
galaxies. In future work we plan to determine the influence of 
other ISM parameters and perform more global simulations of 
these galaxies. 
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